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1. iniiiokk:tion 


The TASC Magsat investigation covers an area in the 
eastern Indian Ocean containing several najor hath^sietric and 
tectonic features (Fig. 1.1-1). The overall objectives of 
this investigation are: 

a Production of siagnetic an<Kaaly maps tr<m Magsat data 
covering the study region (0-50 degrees S, 75 - 115 
degrees E) . 

• Ccmparison of Magsat and satellite altiaeter data 
(Measuring geoid undulation) in this area, and 
quantification of their relationships; 

• Deteraination of the optiaua resolution of Magsat 
ancmaly naps of the study region. 

• Interpretation of Magsat data using satellite 
altiaeter and other geophysical data in order to 
deteraine the origin and sources of the observed 
Magnetic ancauilies. 


The investigation during this quarter focused on 
three areas: 1) an assessaent of the data quality • especially 

characterization of spikes, 2) an analysis of the spectral 
passbands of significance to geophysical analysis, including 
c(»putations of spectral coherence between nearby repeating 
Magsat tracks, and 3) a spectnu Modeling effort designed to 
show the effect of spacecraft altitude on the observability of 
Magnetic an<Hialies. 
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2. SUMMARY OF PROC^ESS 

The work during this quarter eaphasiaed analysis of 
the data quality and a further study of the spectral 
characteristics of the data. During the upc<wing quarters, 
the data will be processed to generate an an<»aly SMp of the 
project area using equivalent source dipoles. This an<Mwly 
nap will then be c<»q>ared to gravity data derived fr<» 
satellite altiaetry. The results of the present iiwrk will 
help in the design of data preprocessing and the choice of 
optiauB grid spacing for ancMaaly aap production. 


• CORRECTIONS FOR DATA GAPS (Section 3.1) 

Data gaps and spikes were identified and corrected. 
Bad data values occuring in the form of known sissing 
data (flagged with a value of 99999.0) were corrected 
with a two'^point interpolation. Snaller spikes 
evident to visual inspection are found to exist after 
rem>val of a field iK>del. 

• SPECTRAL EFFECTS OF REMOVING FIELD MODEL. 

MEAN. AND RAMP (Section 3.1) 

The effect upon spectral estiaates of reaoving a 
field Bodel , aean (D.C. bias), and raap (long 
wavelength c<Msponents) was investigated nuaerically. 
The effects of these preprocessing scheaes were 
analysed for the iapact to analysis of crustal 
anc^aly fields. 

a CCNiERENCE BETWEEN ADJACENT TRACKS (Section 3.2) 

The repeatability of data along tvo satellite passes 
with closely-spaced ground tracks was studied to 
define the wavelength passband containing data of 
geophysical significance. The spectrua of the 
preprocessed data was then analyzed to deteraine the 
coherent bandwidth and noise floor. It was found 
that the spectrua of Magsat along- track data aay be 
divided into 3 aain passbands at wavelengths shorter 
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thiui those roMsved using the field w>del: 1) idiere 

the deta frmi nearly repeating tracks have hig^ 
coherence (i.e. a wavelength band i^ere the data are 
repeatable and considered to be geophysically 
significant, 2800 kat to 800 ka), 2) t^ere the data 
have Xg» coherence but nay have geophysical 
significance (800 lot to 250 ka), and 3) tdiere the 
data has no geophysical significance with regard to 
crustal an(»aly fields ( < 250 ka). 

e ALONG-TSACK SPECTRUM MODELING 
(Section 3.3) 

In order to understand the physical significance of 
along- track Magsat p<Mier spectra, a siaple analytic 
aodel of these spectra was developed. In particular, 
a aodified Bessel function aodel easily provides for 
uiM^ard continuation of along-track Magsat spectra to 
deronstrate their dependence on altitude. 


These analyses yield a better understanding of the 
data and will result in a determination of optimum processing 
techniques to be used in the upc<Hiing work. The r^ainder of 
this report will discuss the details of each of the above 
items. Also during the past quarter, w>rk began on a manuscript 
for sumission to Geophysical Research Letters . Hie results 
given in that manuscript will be more recent and may supersede 
some results presented here . 
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3. TEOPTICAL RESULTS 


3.1 DATA QUALITY AND PREPROCESSING 

Before large scale processing of Magsat data can be 
acco^>lished, it is necessaiiy to evaluate the data quality and 
deteraine what preprocessing, if any, is required. Optiwia 
ccmputational aethods and parameters <e.g. the ainiaua spacing 
for equivalent source dipoles) can then be determined. Iliis 
section will discuss the effects of data preprocessing on spec- 
trua results. 

Twenty- four revs of 80 point averages of scalar aagni- 
tude data crossing the project area (Fig. 1-1.1) were analyzed. 
All revs had 180 points or fewer, and none had fewer than five 
missing points (flagged with 99999.0). In order to maintain a 
uniformly spaced data series, the actual value of these points 
must be estimated and inserted into the time series. For the 
present, two-point linear interpolation was used to estimate 
the missing data values. The time feries were then processed 
to remove a 14th order field model (MGST4/81), and spikes of 
up to 15 gamias are found to contaminate the residual data. 

It is not likely these smaller spikes are real, as their spatial 
w .velengths are shorter than the resolution limit for the data 
(Ref. 1). Since a universally effective algorithm for the 
identification, discrimination and replacement of single-point 
spikes has not yet been developed, these were left in place 
for the time being. Numerical tests have shown that their 
presence does not significantly change the conclusions pre- 
sented here. 
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Effective eaelysie of cruetel fields requires 

a kiKM^ledge of both the wavelength passband i^ich contains the 
desired geophysical information and the effective passband 
where significant information exists in the data. Ehcmledge 
of these passbands will permit determination of the optimum 
spacing for the grid of equivalent dipole sources. As a first 
step in defining these passbands, the effects of r&np r«K>val 
on spectral estimates will be shmm. 

The emmon practice of preprocessing Nagsat an<MMly 
data by removing a mean and ramp alters the spectral content 
of the data. Subtracting a field model of degree and ortler 14 
removes wavelengths longer than about 2800 kil<»eters, primarily 
attributable to the core field. Since this procedure removes 
the core field, the residuals are crustal magnetic an<muilies 
plus external field effects. The external field effects are 
difficult to model, but are at longer wavelengths than crustal 
antmuilies. External field effects can s<»etlmes be removed by 
subtracting a mp fr<» a pass of the residual data. The effect 
of removing a mean or mean and raa^ is illustrated in Figure 
3.1-1. 

Further illustration of mean and ramp reroval is shown 
in Figure 3.1-2. The slope in the midsection of the spectra, 
and the noise floor are unchanged because mean and raaq> r«K>val 
are a form of high-pass filtering. 

Subtracting the mean fr<»i the data rosoves any bias 
fdiich may be due to instrument error or a uniform external 
field. Subtracting a rMq> fr<mi the data is expected to r«K>ve 
primarily the effects of external current fields. Ifowever, 
ramp resK>val will also take out any long-wave length regional 
c<»ponents that may exist in the crustal field. 
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Figure 3.1-1 


Rev 645 8b<H#iiig: e) The eceler Mgnitude date ainua 
the field M>del. b) The autoregressive spectral 
estimate of order 5. The spectra for the atove 
data (curve 1), siean moved (curve 2), and Man 
and rasp renoved (curve 3). Note there is no 
si^ificant difference Mong the 3 spectral 
estiaates at frequencies greater than about 
0.03 Hz (1050 la). Considering the confidence 
bounds of these different estiaatest they are in 
agreMent for f#avelengths longer than about 
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Figure 3.1-2 


Rev 645 is here sh<nm with the Nsgsat date on 
left and Magsat sinus field Model on right. 
(Note the different vertical scales.) Top to 
bott<» rows are: unaltered data, sean r«K»ved, 
and aean and r«sp r«scwed. The dotted lines 
bound the 95% confidence level. 
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The conclueloa it that, for the trecke etomi in Fig. 
1.1-1 , the wavelength paaaband fr<» 2800 ka to 1050 ka 
(extending to ahorter wavelengtha thim contained in the field 
iKHiel) ia edified by the r^Mval of a naan or a Md a 
r«q>. *nie noiae floor appeara to be reached at wavelengtha 
alii^tly ahorter than 250 ka. Thua, the region from 250 ka to 
1050 ka containa inforaation i^ich aay pertain to the ge<^hyaica 
of the legion and ia unalter^ by the roKwal of a aean or a 
Mm md a rMq>. 'Hie tradeoff invoked by tmm or rai^ rmoval 
involvea decreaaing or eliainating the effecta of external fielda 
at the expenae of perturbing the data with wavelengtha Imger 
thm about 1050 ka in the TASC inveatigation area. 


3.2 SPECTRAL ESTlMATlC»f TEOCaQUES iOlD O^f^MCE 

Magaat tracka ifhich are geographically aeparated by 
leaa thm about 100 ka along the entire track aegnent are 
conaidered repeat tracka. Som of theae are ahown in the aap 
of Fig. 1.1-1. Thia separation, md the altitude variation, 
are both mall c<aq»ared to the geophysically ai^ificant 
correlation diatmces in aajor crustal aa^etic moaalies. 
Vithin the usable wavelength paasbmd (1050 to 250 ka) 
described in Section 3.1, c^^tations were perfocmd to find 
idiere the data are repeatable froa track to track (hi^ 
coherence), and «diere they are not (low coherence). 

In thia section, the concept of repeatability ia 
qumtified md illustrated in the frequmcy doaain by 
c<aq>utationa of spectral coherence between nearly repeating 
tracka. Ibeae c<aq>utationa have defined ti#o ba^a of 
differing behavior within the usable data paasbmd (1050 to 
250 ka). A hi^ coherence paasbmd (1050 to 800 ka) defines 
the region where the data are repeatable and the apectrtai of 
the data cm be considered to cmtain geophysically 
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•Ignificaiit iaforMtion. A low coherence peesbead to 250 
la) define* that pert of the date which is not repeatable frcM 
txack to track, bat which m&f still ccmtain aseful data along 
single tracks. The low coherence in this passbi^ is partly 
due to track separaticm «id My also be due to colored <f~*^) 
Mise. 


A quantitative Masure of repeatability is provided 
by the spectral coherence function (Sef. 2) lAich is the 
squared correlation coefficient of the data along the tracks 
at each frequency (or wavelei^th) in the si^ctnni. Spectral 
coherence, p(w), is defined as follows: 


p(w) = 


***<•>*«<■> 


<3.2-1) 


idiere, for the two tracks x(t) and y(t), and are the 
correspcmding auto«spectra and is the cross-spectnai. 
ihatoregressive Mthods of spectrwi estiMtion are discussed in 
Kef. 3. 


'Hie coherence bet%reen pairs of repeat tracks is reduced 
frcui the Mximai value of 1 by non-repeating errors or iw>ise 
in the individual Masurownts. The observed Magsat data is 
coMosed of the siipaal (considered here to be crustal anMMlies) 
and noise (e.g. instnaent noise, orbit errors, tiM-varying 
external fields). The noise is the Mjor cMtributor to the 
iw>n- repeatable portions of the spectra. 

An autoregressive spectral <fStiMti<m algoritlss was 
ai^lied to the sections of Rev 645 (previ<msly described) and 
Rev S61. These ti^ tracks are separated by 49 to 92 km in 
surface distance, and by 70 to ^ loi in altitiwie. Itiis 

1 

I 
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spectral astiaatim algorltte differs froa tliat used ia Sec- 
ti«i 3.1 ia ^at it uses iaforaaticm fi^ both tracks to esti- 
aate the spectra of either track. Thus, i^ere the two tracks 
ate cohereat, the ti«o-ch«^el spectral estiaate is batter (as 
it operates on i^re data). This algoritlai thea estiMtes the 
cc^ereace betaeea the ti#o tracks (eq. 3.2-1), ahich is a aeasure 
of hoa 1^11 oae track cim be predicted fzoa the other. The 
Insult is stoim ia Figure 3.2-1. The aiuaMlp profiles are 
visually colorable, aad they both have a feu spikes (occuriag 
at differeat sai^le poiats). E^aral of the spikes causes the 
estiaated spectra to agree aore closely, but does not sipiifi- 
caatly ch«ige the coaclusioas that aay be drmm irom the ct^ereace. 

'^e spectral structure ia Fig. 3.2-1 is aore detailed 
than shoua «id discussed ia Sectioa 3.1 aaialy because the 
spectral estiaati<m algoritto chose order 10 as the best fit 
here, versus order 5 chosea previously. This aay bet partially 
due to the iacrease of the available data <tuo tracks). A 
spectral peak at about 1700 ka is clearly evident here, md 
uas TOt seen with the order 5 aodel. The cause of this peak 
is a nearly periodic structure of the data, as aay be seen in 
the tiM series of Fig. 3.2-1. 

Choosing a squared coherence level at 0.5 as an arbi- 
traxy cutoff betiMea hi^ coherence and low coherence data, 
the short wavelength liait for coherent data occurs at about 
SOOloi. The long wavelength liait for coherent data (the aaxi- 
MB coherence occurs at about 1^0 ka) is not c<msidered signif- 
ic«it, since unlike the short wavelength liait, it is very 
different Mong different pairs of repeat tracks, and it falls 
within a spectral band affected by roMval of a Meld aodel, 

Man ami ra^». 
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Vork is in progress to add confidence bounds to the 
c<Aerence estinetes. In the future, the ccmfidence bounds 
will be used to define the wavelength band in %^ich there is 
sipiificant coherence. 

The noise floor in these spectra is unchanged fr<m 
that observed in Section 3.1. The data used for these 
analyses i^re taken during Magnetically quiet tines. If the 
na^etic enviromiMiit varied significantly bet%^en the repeat 
tracks, the added noise tmuld be incoherent, thus degrading 
the overall coherence. 

The conclusion fr<» this section is that wavelengths 
longer than alK>ut 800 ka are highly coherent, %diile shorter 
ones are not. Figure 3.2'2 suanarizes these observations 
about the usable inforaation content of each wavelength band. 
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Figure 3.2-2 


Properties of different regions of the Magsat 
along- track pwer spectrua, as deterained frcai 
spectnu analysis and coherence studies. 


The low coherence in the band between 800 and 250 ka 
does not aean that geophysical inforaation cannot be obtained 
in this band. This could be due to variations in altitude and 
geographic position of the repeat tracks. The properties of 
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this band aust b« investigated further with gridded data and 
by ccaqparison with gravity data. 

The results of this and the previous section suggest 
that l<Htf-pass filtering and decimation of Nagsat data into 
saapling intervals substantially longer than the 80 saaple 
averages now available on the Investigator tapes would not 
degrade the analysis of crustal magnetic an<»alies. LcMir-pass 
filtering and decimation over as many as 500 to 600 samples 
could result in more cost-effective data analysis as it «^uld 
only delete incoherent noise. 


3.3 SPECTRUM MODELING 

The along- track power spectrum of Magsat data has 
been estimated and used as a tool to study the data quality 
and resolution capability. It VK>uld also be useful, for 
interpretation purposes, to know idiat these spectra are 
expected to look like for reasonable physical models of the 
earth's magnetic field. Analytic models for these spectra 
could allow study of the effects of instrument noise, external 
fields, and change in altitude. In this section, s<»e 
preliminary results of model develofMient are described. The 
purpose of this «iork was to develop a siJ 4 >le analytic model of 
along- track spectra that could both describe the data and be 
easily upward continued, in order to dei^nstrate the magnitude 
of the effects of satellite altitude on the power spectrum, 
niis is important because features in the spectnmi representing 
crustal magnetic an<»ialies are attenuated rapidly with altitude 
and are obscured by noise when their power drops below a certain 
level. Further work is necessary to relate the parameters of 
this initial model to physical quantities. 
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Upward continuation of a spherical harmonic power 
spectrum (the defined by Lowes, 1974, Ref. 4) is accomplished 
by multiplication of the individual terms of degree n by the 
factor (a/r)^^*^^^^ (where a is the earth’s radius and r is the 
radius of the observation point). For small values of alti- 
tude and for high degree, this spherical harmonic power spectrum, 

R„ = (n+1) I l(g“)^ + (h“)^J (3.3-1) 

" m=0 ^ " 


asymptotically approaches a continuous flat-earth spectrum. 

For cases where the flat-earth approximation is valid (r-a<<a, 

n>>l), a two-dimensional power spectrum, (‘«v upward 

X y 

2 2 1/2 

continued by multiplying by the factor exp[ -h(io +ut„) ' ). 

Thus , ^ y 


‘f'(Wx,Wy,r) = <t>(uij^,u)y)*exp(-2hfl) 


(3.3-2) 


where h=r-a is the height, and fi 



2 

U) +U) 
X 


2 

y* 


For along-track, or great-circle spectra, (such as 
would be observed along tracks of aeromagnetic or satellite 
data), upward continuation is not so straightforward. Formulas 
relating great-circle spectra to spherical harmonic spectra, 
and along-track spectra to two-dimensional isotropic spectra, 
have recently become available but have not been widely used. 
Thus, for example, the along-track (great-circle) aeromagnetic 
spectrum of Alldredge et al . (Ref. 5) was not properly inter- 
preted in many of the older papers which tried to relate these 
observations to global magnetic models. Recently, McLeod and 
Coleman (1980, Ref. 6) derived expressions for the relationships 
between great-circle and spherical harmonic spectra and applied 
these to model satellite observations. 
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To derive a simple possible model of Magsat along- track 
power spectra which permits upward continuation, the following 
assumptions are made: 1) the flat-earth approximation is valid 

for features in the wavelength range of interest, 2) the two- 
dimensional power spectrum is isotropic and has the form of an 
exponential. These assumptions allow the development of an 
analytic expression to describe the along- track power spectrum 
at any altitude. The second assumption is equivalent to stat- 
ing that the two-dimensional power spectrum at any altitude 
may be described by the following function: 

4>(Q) = c exp(-b0)*exp(-2fih) (3.3-3) 


where c and b are parameters, h is altitude, and Q is radial 
frequency for the two-dimensional isotropic spectrum. This 
model is an "attenuated white noise" (AWN) type model since 
<l>(fi)=c at a depth of h=-b/2 (Ref. 7). 


It can be shown that the along-track power spectrum, 
S(u>), corresponding to this two-dimensional power spectrum is 
the Abel Transform of Eq. 3.3-3 (Ref. 8). Thus, 


S(w) 



(3.3-4) 


The solution of this integral for the AWN model is: 


S(u)) = ~ w Kj^[u)(b+2h)] 


(3.3-5) 


where is the modified Bessel function (Ref. 9). 

A choice of values for the parameters b and c was 
made by fitting equation 3.3-5 to an observed Magsat power 
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spectrum at h = 300 km. (The modeling discussed here is being 

done for the spectrum of |B|.) The values obtained were c = 
o ^2 

1.0 X 10 and b = 6.37 x 10 . With these parameter values for 
b and c, the function of eq. 3.3-5 iJ plotted in Fig. 3.3-1 
for altitude values of h=150, 300, and 450 km. If it is assumed 
that the noise floor, representing white noise, is at the 
level of approximately 0.2 (y) /(cy/sample), then the resolution 
limit of satellite data would be approximately 160 km for 
h = 150, 350 km for h = 300 km, and 540 km for h = 500 km. 


SPECTRUM OF |B| ANOMALY, ORBIT 1201 



3600 1200 360 72 


WAVELENGTH (KM| 


Figure 3.3-1 Observed Magsat spectrum for |B| of orbit 1201 

compared with values given by the analytic 
model of Eq. 3.3-5. The effect of different 
heights h, is shown. (Orbit 1201 was at h=350 
to 360 km during the pass over the investigation 
region. ) 
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This model has certain simplifications: most notably 

the flat-earth approximation and isotropy of the two-dimensional 
power spectrum. However, the estimates of the resolution capa- 
bility at different altitudes seem reasonable. The formula of 
eq. 3.3-4 is very similar to the formula of eq. 27 of McLeod 
and Coleman (1980, Ref. 6), which is the asjmptotic approximation 
of an expression describing the relationship between an average 
great-circle power spectrum and the spherical harmonic spectrum. 
Further work will give the as 3 nnptotic form of the spherical 
harmonic spectrum, R^, which is implied by the flat-earth model 
of eq. 3.3-5. 


When plotted on a log-linear scale, the function (eq. 
3.3-5) has the form of a straight line with a slope=-2h, indi- 
cating that it is essentially exponential. This is the foim 
of spectra observed in aeromagnetic surveys as modeled by Spec- 
tor and Grant (1970, Ref. 10) for the purpose of determining 
depth to magnetic basement. The fact that our along-track 
spectrum estimates from Magsat data show power law, rather 
than exponential behavior, needs further investigation. This 
difference is important because many models have assumed the 
exponential form. Reasons for this disagreement could be: 1) 
the crustal field is not exponential and 2) time-varying noise 
sources . 


PLANS FOR THE NEXT REPORTING PERIOD 

The Magsat data will be prepared for computing crustal 
anomaly maps by modeling the equivalent dipoles in the project 
area. Special emphasis will be placed on the removal of spikes 
and proper filtering and preprocessing of the data. Since two 
more investigator tapes have been recently received (Sept. 

30), production processing can begin on the complete data set. 
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The first runs of the equivalent source fitting program will 
be based on a 7 degree spacing (as data is highly coherent to 
a short wavelength limit of 800 km - Section 3.3). 


COST DATA 

Total expenditures through 30 September have been 
$23,426. Expenditures during the quarter 1 July through 30 
September were $12,093. 
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